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The Laplacian A (^^ + in the plane) is one 
of the most basic operators in all of mathematical 
analysis. It can be used to construct the impor- 
tant spacetime equations of mathematical physics, 
O ^uch as the heat equation, the wave equation and 
D the Schrodinger equation of Quantum Mechanics. 

It has been studied from many points of view, and 
'xj" in many different contexts (Riemannian manifolds, 

graphs, fractals, etc). The eigenvalue equation 
p ^ —Am = Am in a domain with suitable bound- 
ary conditions (Dirichlet conditions — and 
. Neumann conditions f^laa = are the most com- 
1-Q mon) defines both the eigenvalue A (in many cases 
d a nonnegative real number) and the eigenfunction 
^ u. The set of all eigenvalues (usually a discrete 
I— iset) is called the spectrum of the Laplacian on fl, 
and there is a vast literature on the relationship 
between the spectrum and the geometry of fl. See 
^ for example jK], |GWW] and [Z2]. The spectrum 
00 contains a lot of information about the domain, 
r — and teasing out this information, including rela- 
tionships between quantum and classical mechan- 
^\ ics on the domain, has been a fascinating, ongoing 
and highly nontrivial enterprise. For example, the 
famous Weyl law gives the asymptotic size of the 
. . kth eigenvalue as fc — > oo. There are also inter- 
^ esting and intricate estimates on the size of the 
first and second nontrivial eigenvalues. Needless 
to say, there are no connections between estimates 
of small and large eigenvalues. 

Eigenfunctions are also fascinating objects. 
Clearly these are more complicated objects, so it 
can be expected that it is more difficult to say 
things about them. Many of the classical special 
functions and orthogonal polynomials can be re- 
lated to eigenfunctions of the Laplacian for specific 
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geometries. There have been many recent stud- 
ies concerning eigenfunctions associated to large 
eigenvalues. In this note we invite you to look at 
some startling pictures of some specific "localized" 
eigenfunctions associated to small eigenvalues. 



Examples of Localization 

How should we define a "localized" eigenfunc- 
tion? We would be tempted to say it is an eigen- 
function with support ili that is considerably smaller 
than all of il. But it is well-known that eigen- 
functions are real analytic functions, hence cannot 
vanish on any open set. So we must be content 
with saying that the function is "very small" on 
the complement of fii. This is of course not a 
mathematical definition, although it might be ac- 
ceptable to a physicist, or a Justice of the Supreme 
Court. One could make it into a precise definition 
with a parameter e to quantify the statement "very 
small" , but this just begs the question: how small 
does e have to be to make the statement interest- 
ing? In this note we will show you some pictures 
to try to convince you that there are surprising ex- 
amples where e is smaller than you might expect. 

Localized eigenfunctions have been observed 
before. As usual, physicists know more than math- 
ematicians in the subject, but with less certainty 
[this is a kind of uncertainty principle]. Regardless, 
the cross-pollination in this subject between these 
two groups over the past century merits recognition 
and esteem. For high frequency eigenfunctions, re- 
lations between eigenfunction localization and bil- 
liard dynamics have been studied. A nonexhaus- 
tive list includes treatments of (non)localization 
on: closed stable geodesies |BLj and closed unsta- 
ble geodesies [HI [CH IH [HH]. Other results ad- 
dress dichotomies [Be], numerical aspects |Bal IBaj . 
rarity |Sn[ ICll IZlj and near-nonexistence [L] of 
such (phase space) localized eigenfunctions, as A — > 
oo. In the low frequency realm, no deep expla- 
nation for eigenfunction localization seems to ex- 
ist. Low frequency, or "ground state" eigenfunc- 
tions have been widely studied (see for example 
[P] IBMP| . However, the authors can only find 
scattered examples of low frequency localization, 
such as: near a fractal boundary |RSH] . in narrow 
channels between domains |CHj and in square pairs 
with irrational ratios of frequency oscillations (Ex- 
ample 3 in |JMSj ). As is well known, an eigenfunc- 
tion is itself an eigenfunction in each of its nodal 
domains (with appropriate boundary conditions). 
Therefore, results for low frequency eigenfunctions 
can inform higher frequency studies. 
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Figure 2. Two rooms with a passage 
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Figure 1. Localized Circle Eigenfunction 

A simple example of an eigenfunction on the 
disc localized to a neighborhood of the boundary 
circle is shown in Figure [T] These examples tend 
to involve eigenfunctions with eigenvalues rather 
high up in the spectrum and domains with spe- 
cial types of billiard flows. In contrast, our exam- 
ples occur low down in the spectrum and are con- 
sequences of symmetry considerations. We work 
with Neumann boundary conditions because they 
are natural (the weak formulation of the eigen- 
value equation is — /^(Vm- Vw) dx = X J^^ uv dx for 
all reasonable test functions v, without imposing 
any boundary conditions) and they were essential 
in our work on approximating fractal Laplacians 
with ordinary planar Laplacians [BHS], |HSj . Sim- 
ilar examples with Dirichlet boundary conditions 
also exist. It was our coauthor Tyrus Berry who 
first serendipitously discovered examples of local- 
ized eigenfunctions on sawtooth shaped domains 
as reported in |BHS| . but these examples did not 
play any role in the theory developed there. By 
coincidence, many localized eigenfunctions of frac- 
tal Laplacians have been known since the work of 
Fukushima and Shima |FS] . and these can also be 
explained by symmetry considerations |BK| . See 
|Stll ISt2] for expository accounts of this phenom- 
enon. 

Our examples can be though of as modified 
versions of "rooms and passages" domains |CHj . If 
il = fii U U consists of two rooms, fii, ^2, 
and a very short narrow connection passage U,^ (see 
Figure |2| then it is not surprising that there are 
Dirichlet eigenfunctions on VI that are very close 
to Dirichlet eigenfunctions on fii extended to be 
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Figure 4. Cow Domain, height 
h = 0.1 

zero on U fis. In our examples we will join 
two rooms = f2i U ri2 that intersect in a small, 
but not very small piece. Figure [3] shows one ex- 
ample, the "smiley face", and Figure |4] shows an- 
other, the "cow". The main idea is that fii must 
possess an axis of symmetry, and the rooms join 
together at one end of this axis. In other words, 
there is a line L such that the reflection R in L 
preserves fii, i?f2i = fii, and L passes through 
r^i n ri2- Suppose ui is a Neumann eigenfunc- 
tion of the Laplacian on fii that is skew-symmetric 
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with respect to R, ui{Rx) — —ui{x). Such eigen- 
functions occur throughout the spectrum, since in- 
deed every eigenspace spHts into functions that 
have symmetry and skew-symmetry with respect 
to R (most eigenspaces are 1-dimensional and are 
of one or the other symmetry type). Then ui van- 
ishes along L. If the point p where L intersects 
dCli is a corner point, then ui and Vui vanish at 
p, so Ui is relatively small near p. (For example, 
ui{x, y) = cos TTX — cos ■ny at the origin if ili is the 
unit square.) So it is not surprising that there is 
an eigenfunction m on 17 that is close to ui on f2i 
and close to zero off 51 1. 

Such reasoning does not yield a sharp estimate 
for how localized w is, so we look at some exper- 
imental evidence. We use Matlab to numerically 
approximate some eigenfunctions on our domains 
using the finite element method. Figures [5]{T2| show 
the results. Note that we normalize the eigenfunc- 
tions to have norm on f2 equal to 1, and we can 
measure the localization either by the norm on 
\ ill, a kind of average localization, or the L°° 
norm on \ f2i, a uniform localization. In Figures 
|6|7|9| and [lO] we show both of these as a function 
of the aperture size [height., will edit later] for the 
connection for the each domain. These are log-log 
plots, suggesting a power law relationship over the 
given range of h values. In the tables below we 
give the best fit power law for two eigenfunctions 
in each domain. Note that the powers vary consid- 
erably in these four examples. 
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Conclusion 

Our examples show how easy it is to find sur- 
prisingly localized eigenfunctions. The domains do 
not have to have any special properties beyond the 
symmetry of the f2i piece (breaking the symme- 
try even slightly makes the examples disappear). 
We do not have to go very high up in the spec- 
trum. As mathematicians it is natural for us to 
want a theorem that explains the examples, or at 
least a conjectured theorem. Perhaps there is such 
a theorem, and a perceptive reader might be able 
to find one, but at present we don't see any. Or, 
we might suggest that there is more to mathemat- 
ics than just theorems. This might sound like a 
radical suggestion, or perhaps it is just common 
sense. 




Figure 5. Smiley, Fifth Eigenfunction 




Figure 6. Localization in 
norm, for Figure |5] example 




Figure 7. Locahzation in L' 
norm, for Figure [5] example 
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Figure 8. Cow, Fourth Eigcnfunction 




Figure 9. Localization in 
norm, for Figure [9] example 






Figure 11. Smiley, Twelfth Eigenfunction 
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Figure 10. Locahzation in L°° 
norm, for Figure [9] example 



Figure 12. Cow, Eleventh Eigenfunction 



LOCALIZED EIGENFUNCTIONS: HERE YOU SEE THEM, THERE YOU DON'T 



5 



References 

[BL] V. M. Babic and V. F. Lazutkin, The eigenf unctions 
which are concentrated near a closed geodesic. 1967 
Problems of Mathematical Physics, No. 2, Spectral 
Theory, Diffraction Problems (Russian), p. 15-25 Iz- 
dat. Leningrad. Univ., Leningrad 

[Ba] A. Backer, Numerical aspects of eigenvalue and eigen- 
function computations for chaotic quantum systems. 
The mathematical aspects of quantum maps, 91-144, 
Lecture Notes in Phys., 618, Springer, Berlin, 2003. 
MR2159324 

[BMP] R. Banuelos and P. Mendez-Hernandez, Sharp in- 
equalities for heat kernels of Schrodinger operators 
and applications to spectral gaps. J. Funct. Anal. 176 
(2000), no. 2, 368-399. MR1784420 (2001f;35096) 

[Ba] A. H. Barnett, Asymptotic rate of quantum ergodic- 
ity in chaotic Euclidean billiards. Comm. Pure Appl. 
Math. 59 (2006), no. 10, 1457-1488. MR2248896 
(2007d:81085) 

[BK] M. T. Barlow and J. Kigami, Localized eigenf unctions 
of the Laplacian on p.c.f. self-similar sets. J. London 
Math. Soc, 56 (1997), no. 2, 320-332. MR1489140 
(99b:35162) 

[Bo] M. V. Berry, Regular and irregular semiclassical wave- 
functions. J. Phys. A 10 (1977), no. 12, 2083-2091. 
MR0489542 (58 #8961) 

[BHS] T. Berry, S. Heilman, and R. S. Strichartz, Outer 
Approximation of the Spectrum of a Fractal Laplacian. 
(to appear) Exp. Math. arXiv:0904.3757vl 

[BZ] N. Burq and M. Zworski, Bouncing ball modes and 
quantum chaos. SIAM Rev. 47 (2005), no. 1, p. 43-49 
(electronic) . 

[CI] Y. Colin de Verdiere, Ergodicite et fonctions pro- 
pres du laplacien. [Ergodicity and eigenfunctions of the 
Laplacian], Comm. Math. Phys. 102 (1985), no. 3, p. 
497-502. MR0818831 (87d:58145) 

[C2] Colin de Verdiere, Yvcs(F-GREN-F); Parisse, 
Bernard(F-GREN-F) Equilibre instable en regime 
semi-classique. L Concentration microlocale. [Semi- 
classical unstable equilibrium. I. Microlocal concentra- 
tion] Comm. Partial Differential Equations 19 (1994), 
no. 9-10, 1535-1563. MR1294470 (96b:58112) 

[CH] R. Courant and D. Hilbert, Methods of mathemati- 
cal physics, vol. IL Interscience Publishers, New York, 
1962. MR0140802 (25 #4216) 

[HH] A. Hassell and Luc Hillairet, Ergodic billiards 
that are not quantum unique ergodic. (preprint) 
arXiv:0807.0666v3 

[HS] S. M. Heilman and R. S. Strichartz Homotopies 
of Eigenfunctions and the Spectrum of the Lapla- 
cian on the Sierpinski Carpet, (to appear) Fractals. 
larXiv: 0908.2942vl 

[H] E. J. Heller, Bound-state eigenfunctions of classically 
chaotic Hamiltonian systems: scars of periodic or- 
bits. Phys. Rev. Lett. 53 (1984), no. 16, 1515-1518. 
MR0762412 (85k:81055) 

[FS] M. Fukushima and T. Shima, On a spectral analysis 
for the Sierpinski gasket. Potential Anal. 1 (1992), 1- 
35. MR1245223 (95b:31009) 

[GWW] C. Gordon, D. L. Webb and S. Wolpert, One 
cannot hear the shape of a drum. Bull. Amer. Math. 
Soc. (N.S.) 27 (1992), no. 1, p. 134-138. MR1136137 
(92j:58111) 

[JMS] P. W. Jones, M. Maggioni and Raanan Schul, Man- 
ifold parametrizations by eigenfunctions of the Lapla- 
cian and heat kernels. Proc. Natl. Acad. Sci. USA 105 
(2008), no. 6, p. 1803-1808. 



[K] M. Kac, Can one hear the shape of a drum? Amer. 

Math. Monthly 73 1966 no. 4, part II, p. 1-23. 

MR0201237 (34 #1121) 
[L] E. Lindenstrauss, Invariant measures and arithmetic 

quantum unique ergodicity. Ann. of Math. (2) 163 

(2006), no. 1, p. 165-219. 
[P] Payne, Lawrence E. On two conjectures in the fixed 

membrane eigenvalue problem. Z. Angew. Math. Phys. 

24 (1973), 721-729. MR0333487 (48 #11812) 
[RSH] S. Russ, B. Sapoval and O. Haeberle Irregular and 

fractal resonators with Neumann boundary conditions: 

Density of states and localization. Phys. Rev. E 55 

(1997), no. 2, p. 1413-1421. 
[Sn] A. I. Snirel'man, Ergodic properties of eigenfunctions. 

Uspehi Mat. Nauk 29 (1974), no. 6(180), p. 181-182. 

MR0402834 (53 #6648) 
[Stl] R. S. Strichartz, Analysis on fractals. Notices Amer. 

Math. Soc. 46 (1999), no. 10, 1199-1208. MR1715511 

(20001:58035) 

[St2] R. S. Strichartz, Differential Equations on Fractals: 
A Tutorial. Princeton University Press, Princeton, NJ 
2006. MR2246975 (20071:35003) 

[Zl] S. Zelditch, Uniform distribution of eigenfunctions on 
compact hyperbolic surfaces. Duke Math. J. 55 (1987), 
no. 4, p. 919-941. MR0916129 (89d:58129) 

[Z2] S. Zelditch Inverse spectral problem for analytic do- 
mains, II: Z2-symmetric domains. Ann. of Math. (2) 
170 (2009), no. 1, 205-269. 



